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Abstract 

We study the late-time behaviour of tilted perfect fluid Bianchi type 111 models using 
a dynamical systems approach. We consider models with dust, and perfect fluids stiff er 
than dust, and eludicate the late-time behaviour by studying the centre manifold which 
dominates the behaviour of the model at late times. In the dust case, this centre manifold 
is 3-dimensional and can be considered as a double bifurcation as the 2 parameters (h 
and 7) of the type VI^ model are varied. We calculate the decay rates and show that for 
dust or stiffer the models approach a vacuum spacetime, however, it does so rather slowly: 
p/H 2 ~l/lnt. 



1 Introduction 

In previous papers the dynamical behaviour of tilted Bianchi models have been studied [TJ GH 
Sl[5] using the dynamical systems approach [6]. In particular, the tilted Bianchi models of type 
VI/j was studied in [4j, where it was pointed out that the type III=VI_i model is a bifurcation 
value in terms of the group parameter h. The analysis of this model therefore requires a more 
complicated centre manifold analysis. In addition, in terms of the equation of state parameter 
7 there is also a bifurcation at 7 = 1, which results in a 3-dimensional centre manifold in this 
casfl The dust type III model is therefore a particularly interesting model as it appears as a 
double bifurcation (this was also pointed out in [5]). Due to the centre manifolds arising from 
these bifurcations, the Bianchi type III models with 1 < 7 < 2 was not studied in detail in [4j. 

The type III model also have an interesting geometric interpretation. The type III Lie 
group can be considered as the Thurston geometry H 2 x R which plays an important role in 
3-dimensional geometry [7] . Interestingly, this Thurston geometry also permits for a type VIII 
action which allows for a possible connection between the behaviours of these models. In [4] 



*Address after April 15, 2008: Dept. of Mathematics and Natural Sciences, University of Stavanger, N-4036 
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1 The 3rd dimension of this centre manifold is related to the existence of a line of equilibrium points and is 
present even for non-tilted models. 
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such a connection was pointed out; however, the oscillatory behaviour present in the type VIII 
models [3] have no analogue in the type III models. 

The tilted Bianchi type III model is the last of the ever-expanding Bianchi models left to 
study in terms of its late-time behaviour. It is the aim of this paper to complete this analysis. 
Some partial results are known; for example, the case < 7 < 1 was analysed in [4] (including 
an analysis in some particular subspaces). Here, we will consider the fully tilted type III model 
and show that for 1 < 7 < 2 the late-time asympotote is the self-similar vacuum spacetime 
given by 

ds 2 = -dt 2 + t 2 (dx 2 + e- 2x dy 2 ) + dz 2 . 

The tilt, on the other hand, depends on 7 in the following way: for 7=1, the tilt tends to zero 
(albeit slowly), while for 1 < 7 < 2 the tilt is asymptotically extreme. 

The dynamical system considered in this paper is a prime example of a dynamical system 
appearing in other theories of applied mathematics and physics. This dynamical system, which 
is also constained, is of relatively high dimension, namely 7. The analysis performed here is 
a certain aspect of this dynamical system, namely a centre manifold analysis. This centre 
manifold appears as two of the parameters of the theory experience bifurcation values. We 
apply centre manifold theory, and for all the variables of the theory, we find the decay rates 
and the correction terms, as the dynamical time t — > 00. 

2 Equations of motion 

The dimensionless variables of the theory are the shear variables, (S+, £12, £13, £23); g eo ~ 
metric variables, {A, N, A); and fluid variables (Q, v\,V2, ^3). The variable f2 is the (expansion- 
normalised) energy-density while Vi are so-called tilt variables (they correspond to the velocity 
of the fluid). If < V < 1, where V 2 = v\ + v 2 + v 2 , we call the model tilted; the special 
cases V = and V = 1 are called non-tilted and extremely tilted, respectively. Moreover, the 
dimensionless time-variable, r, can be related to proper time, t, by 

dT -H 



where H is the Hubble scalar. 

The papers [2J 0] contain all the details regarding the determination of the evolution equa- 
tions for the models under consideration. 
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The equations of motion for the general tilted type III model are: 





= (q- 2)E+ + 3(E 2 2 + E 2 3 ) -2N 2 + ^- {-2v\ + v\ + v\) 


(1) 


E'_ 


= (q-2- 2V3E 23 A)S_ + V3(S? 2 - E? 3 ) + 2A7V + ^f - ( V 2 - W 2 ) 


(2) 


y' 


= (q - 2 - 3E+ - V3£_) E 12 - V3 (E 23 + E_A) E 13 + ^^1% 


(3) 


^13 


= (q - 2 - 3E+ + V3S_) E 13 - V3 (E 23 - E_A) E 12 + ^^ Vl v 3 


(4) 


y' 

^23 


= (g - 2)E 23 - 2\/lW 2 A + 2V3AE 2 . + 2V3E 12 E 13 + ^p^v 2 v 3 

G+ 


(5) 


N' 


= (q + 2E+ + 2V3E 23 A) N 


(6) 


A' 


= 2V3E 23 (1-A 2 ) 


(7) 


A' 


= (<? + 2E + )A 


(8) 



The equations for the fluid are: 

n' = ^-{2q-(3 1 -2) + 2 1 Av 1 + [2q(~f-l)-(2-~ f )- 1 S}V 2 } (9) 

v[ = (T + 2y + ) Vl - 2\/3E 13 v 3 - 2\/3E 12 i; 2 -A(v 2 + vj) - V3N (v% - v 2 ) (10) 

v' 2 = (T-E+-v^E-)^-V^(S23 + S_A)«3 + >/3AJVui«3+(i4 + >/3Jv)wi«2(ll) 

v' 3 = (T-E + +V3E_)w 3 -V3(E 23 -E_A)i; 2 -V3AiVwiw 2 + (A-V3iv)wiw 3 (12) 

v> = i [(37 ~ 4) ~ 2(7 ~ ~ 51 ' (i3) 

where 



3 = 


0^2 , l(3 7 -2) + (2-7)V 2 
2L + 2 l + ( 7 -l)V 2 52 


E 2 = 


E 2 ^ + E_ + E 2 2 + E 2 3 + E 23 


5 = 


E Qb c a c fc , c a c a = 1, w a = w\ 


V 2 = 


t; 2 +u| 


T = 


[(3 7 - 4) - 2( 7 - 1)M] (1 - V 2 ) + (2 - 7 )V 2 S 


1 - (7- 1)V 2 


G+ = 


l + ( 7 -l)y 2 . 
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These variables are subject to the constraints 



= E + A + N +CI 



2E4- A + 2E-N 



Z 12 (N + V3A) + £ i3 AA 

Eia(JV- V3A) + £ 12 AA 
A 2 -3(1- A 2 ) A 2 . 



7^2 

G + 

7^3 



(14) 
(15) 

(16) 

(17) 
(18) 



In this paper we will be concerned with the case when the parameter 7 obeys 1 < 7 < 2. 
The case < 7 < 1 is considered in [4]. The generalized Friedmann equation ([14"]) yields 
an expression which effectively defines the energy density il. We will assume that this en- 
ergy density is non-negative: > 0. Therefore, the state vector can thus be considered 
X = [£+, £_, £12, £13, £23, A, A, A, vi, V2, V3} modulo the constraint equations (fl~5|) - (|T5]) . The 
evolution takes place on a seven dimensional submanifold of the twelve dimensional space; thus 
the dimension of the physical state space is seven. 
Additional details are presented in [2j [4] . 

2.1 Fluid Vorticity 

In general, the fluid will have a non-zero vorticity, W a . The vorticity of the fluid for the type 
III model is given by: 



Wa = ( Nabyb + £abcvbAC + 7-^y2 N bcv b v c v a j . ir„ = -r"ir„. 



(19) 



where 



B 



1 - \{V 2 + V 2 S + 2A a v a ) 
Vl-^ 2 [l-(7-l)^ 2 ] 



2.2 Equilibrium points 

The equilibrium points of the Bianchi type III model with 1 < 7 < 2 are all given in [1]. 
We will be concerned with the general case for which the following equilibrium points are of 
importance. 

2.2.1 B(I): Equilibrium points of Bianchi type I 

1. 1(1): £+ = £_ = £12 = £13 = £23 = A = N = V = and Q, = 1. Here, |A| < 1 and is 
an unphysical parameter. 

2.2.2 T(III): Vacuum case (fi = 0) 

All of these equilibrium points are plane wave solutions and have 

Q = £ 12 = £13 = £23 = 0, £_ = = v / -£+(l + £+), 
A = (1 + £+), - 1< £+ < 0, |A| < 1. 
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It is also advantageous to introduce r = \f\ — A 2 . This implies that we can write 



We will also define p by 



p = v\+v\. 



The equilibrium points are then determined by the tilt velocities: 

1. £(—1): v\ = V2 = V3 = 0. 

o 7l 1\ 7(l+3r 2 )-2(l+2r 2 ) „ 2(l+3r 2 ) ^ ^ „ 

2. £(-1): vt = 2ra( -y-i) S «2 = «3 = 0, \ +5r2 < 7 < 2. 

3. C±(-l): wi=±l,W2=«s = 0. 

4- «* = - 7(1+ £j^ 5r2 \ fe 2 - vl) = pr, 

7 )[ 7 (l+3r 2 )-(l+5r 2 )] 1+5r 
4r 4 (2- 7 ) ' wnele i +3r 

^ < 7 < for r(r + V^TT) < 1. 



P = (1+r2)(3 ' 27 4 ) Il ( 2 1 -:f ) ' (1+5r2)1 - ^ere < 7 < 1 for 1 < r(r + V^TT) and 



5. T (—1): This case collapses to the line-bifurcation 7 = 1, v\ = 0, (w 2 — w 2 ) = — pr, 
< p < 1. 



6. £+(-1): Here, 1 > r(r + Vr 2 + 1) and «i = --p 21 ^, (v| - vfj = pr, p = 1 - v\ = 



(l+r 2 )(l-3r 2 ) 
(1-r 2 ) 2 ' 



7. Here, -1 < r(r - Vr 2 + 1) and ui = 0, (w| - w|) = -r, p = 1. 

2.3 Equilibrium points of special importance 

The following equilibrium points all have non-negative eigenvalues and are therefore potential 
late-time attractors: 



1. V(III): This is the special case T (—1] 
vacuum solutions for h = — 1, 7 = 1: 



and is a special line-bifurcation of 'tilted' 

!-=! 



+4' 4 ' 4' 

£12 = £13 = £23 = A = 0, v\ = v 2 = 0, < v 3 < 1. 



The extreme limit of this equilibrium point, lim„ 3 ^i V(III) = £ p (— 1] 



. We will call 

r=l 

£"(-1) for Sp(IH) for simplicity. We will also define P {III) = \im V3 ^ V(III) = 

A-l)lr=T 

Regarding the eigenvalues; there are 3 eigenvalues which are zero and the rest all have a 
negative real part. In order to resolve the stability properties of these equilibrium points 
one has to resort to centre manifold theory. 

2. P (j/j) = r(-i)| p=1)7 = i. 



3. e-{iii) = e-{-\) 



, 1< 7 <2. 

r=l 
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3 The late-time behaviour 

In [4] the eigenvalues of all the equilibrium points were computed. This analysis showed that 
the only potential equilibrium points acting as local attractor s are V{IIT) and E~{III). Fur- 
thermore, an extensive numerical analysis indicates there are no other kind of attractors (like 
attracting curves). Lastly, the existence of monotone functions [1] also restricts the attractors 
to lie in certain subsets. Although, we have no proof, this indicates that the following local 
attractors are indeed global attractors (except for a set of measure zero) . 

For the 7=1 line of equilibria V{III), two of the zero eigenvalues of the linearised system 
actually correspond to a non-trivial Jordan block; i.e., one of the Jordan blocks is of the form 




1 
" 



This means that a generic solution 'drifts' along the line of equilibria; the amount of 'drifting' 
depends on the second order terms. However, it can be shown that the line V(III) is indeed 
unstable; the solutions drift towards the end point Vo(III). 



3.1 The case 7 = 1: Attractor is V (III). 



The centre manifold in this case is 3-dimensional. The potential attracting equilibrium point 
is Vo(IH)- The centre manifold can be found as follows. It is useful to define £+ and £_ by 
rotating £ + and £_ according to: 



1. 



V3g V3g 
2 ' 2 H 



1. 



-£_ 



(20) 



Let us define the variables x 1 : 

(^•+^23,N,X,V 1 ,V 2 ,V 3 ) 



1 v/3 

- + x 1 ,x 2: — +X3,X4,x 5 ,a; 6 ,x 7 



To determine the late-time behaviour it is actually necessary to expand the variables to 3rd 
order in x = (x l ) (2nd order is not sufficent to determine all the coefficients of the leading order 
terms of the variables). The variables, £_, £12 and £13 can be found by solving the linear 
constraints, eqs. (| 1 5ti 1 T[) . while O can be found by solving the Friedmann constraint, eq. (fT4|) . 
The equations can now be written: 



x' = Lx+ C(x,x) +T(x,x,x) + C(x 4 ), 



(21) 



where C and T is a bilinear and trilinear vector-valued function, respectively. 

Now, it is advantageous to define y = Px so that the linear term in y is the Jordan canonical 
form of L. This can be accomplished by defining 



1 
1 

-1 
-1 












1 
1 

-1 
-1 





2V3 
-2V3 
-2V3 
2V3 










4 



4 






(22) 



THE TYPE III BIFURCATION 



7 



In terms of y we get 

y' = Jy + C(y,y)+T(y,y,y)+0(y 4 ), (23) 

where 

J = diag(o, 0,4 0,-| ,-§ 

From this we note that the 3-dimensional centre manifold can be parameterised by the variables 
{x,y,z) = (2/1,2/1+2/3,2/3 + 2/5)- 

The next step is to find functions Yi(x,y, z), i = 2,4,6,7 so that yi — Yi(x,y,z) — are 
invariant submanifolds. Since the centre manifold will dominate the late-time behaviour, the 
variables 2/i, i — 2, 4, 6, 7 will behave as 2/i ~ Yi(x, y, z) at late times. 

To third order in (x, y, z) we obtain: 

Y 2 (x,y,z) w -y 2 + x{x + z) 2 , 
Y 4 (x,y,z) w y 2 + -yx - x(x + zf , 

Y 6 (x,y,z) w ^x[-j/ + 2(x + z) 2 ], 

y 7 (x,y,z) w -^(x + z)[4x(x + z) + V3y]. 

Substituting this into the equations for (x,y,z) we finally get the equations on the centre 
manifold: 

x' = x 2 + \x z + IQxy 2 + G(y 4 ), 
2/' = 2yx-^yx 2 - 6y 3 + G(y 4 ), 

z' = xz - ix 3 - I3xy 2 -3y 2 z + (D(y 4 ). (24) 
These equations can now be solved to give: 



T - To 3(T-T ) Z 



C 2Cln(r-r ) ,^ r/ , , n 

(t-t ) 2 3(T-T ) d 



T — To 6{T — To)^ 

where r , vq and C are constants. Substituting the solutions into the original variables, we get 
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the decay rates: 



8- - ^M^( 1 + |^) +0( ^,, 

E12 . _^( 1+ ^) +0(T - 5) , 

Sl , . ^( 1+ ^) +0(T -3), 



T 



„ . i( 1+ i^i) +0(O , 

* - + (25) 

The constant, To, has been eliminated by a shift of time, r. We have also included the second 
term in the expansion which is decaying as (In r)r _1 compared to the 1st order term, and the 
order of the correction terms for all the variables. 

The Hubble scalar and the cosmological time can also be calculated to give 

H « H rh-i T , t-to« ^-r-4 e l T . (26) 

This shows that the attractor is indeed the equilibrium point Vq(III). It is therefore 
plausible that all solutions, except a set of measure zero, approach this equilibrium point as a 
power law at late times; in particular, the expansion-normalised energy density goes as f2 cx 1/r 
for sufncently large r. 

Let us consider some other physical consequences of this result. As we now how the decay 
rates, we can calculate the behaviour of certain physical quantities. Consider, for example, the 
fluid vorticity given in eg. ([TO)) . The fluid in this case is pressureless matter (dust). At late 
times, the (square of) the spatial vorticity of a pressureless matter in a Bianchi type III model, 
will decay as: 

W a W a <x-^. 

Hence, the vorticity decays fairly slowly. Therefore, due to the presence of the centre manifold, 
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certain physical quantities will decay slower than the usual exponential decay seen in most 
other Bianchi modelfl 



3.2 The case 1 < 7 < 2: Attractor is £'(111). 

For the case 1 < 7 < 2 the potential attractor is the extremely tilted equilibrium point £~ (III) ■ 
In this case the centre manifold is 2-dimensional so the analysis simplifies slightly. We can 
simplify even further by restricting to the extremely tilted invariant subspace V = 1 in which 
the centre manifold lies. 

The analysis is analogous to the 7=1 case; however, we can also solve for V3 = yl — u| — lif . 
We again define the variables x 1 as follows: 



(Y, + ,Y<23,N, A, v 1 ,v 2 ) 



\ V3 
- +Xl,X 2 ,— + x 3 ,x 4 ,x 5 ,x 6 



We will now expand the equations of motion to 2nd order in x = (x l ). The variables, S12 
and £13 can be found by solving the linear constraints, eqs. (|I5H17[) . while f2 can be found by 
solving the Friedmann constraint, eq. (fM]) . as before. 



We again align our variables with the Jordan canonical form, so we define: 



P = 



1 
1 

-1 
-1 










4V3 
3 



9 

, 8\/3 
T 9 



2\/3 

2V3 

-2^3 

I - 2V3 
2 

_ 5 
3 







f 



_ 2 
3 

2 
3 









_ 4 
3 

4 

3 ■ 



(27) 



The Jordan matrix is 



J = diag 



3 3 3 
2'~2 '~2 



hence, we can use y\ and j/3 to parameterise our centre manifold. 

A similar procedure as in the dust case gives us the following equation of motion on the 
centre manifold: 



Vi 

y's 



8y? + 0(y 3 ), 
&yiy 3 + 0(y 3 ). 



(28) 



To lowest order, this can be integrated to yield y\ = — l/(8r), y 3 — (vq + 1)/(8t), where 
a constant has been eliminated by a shift of time. Substituting these expressions into the 



2 There might be other physical consequences as well, for example in the cosmic microwave background. 
Such a slow decay of the vorticity will most likely have observational consequences for the cosmic microwave 
background radiation |8|. 
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variables gives the decay rates (the errors are also indicated) : 



Sis = -^+C[(lnr)/r 2 ], 

A = |i + G[(ln r)/r 2 ] 

fi = 1 + G[(ln r)/r 2 ], 

«i = + 0[(ln t)/t% 

St 

V2 = -^-+0[(lnr)/r% 
lor 

1 9w 2 + 256 _ r/1 , . 3l 

, / , \ (2-y- 3) 3(7-1) 

a/1-F 2 = (Vl-^ 2 J T -4<2^) e -T2^r^. (29) 

We note that also in this case the energy-density approaches zero, but it does so rather slowly: 
f2 oc 1/t. The tilt, on the other hand, evolves exponentially towards extreme tilt. 
The Hubble scalar and the cosmological time are found to be 

H^H e-i T ri, t - to » ^U* t t"*. (30) 

From this result we can see that the cosmological time, which corresponds to the proper 
time of observers following the geometrically defined congruence, is unbounded as t — > 00. If 
we instead consider the congruence defined by the four- velocity of the fluid [9] , the proper time 
is given by 

f T 1 / 3(37-4) _ 

T-T = —Vl-V 2 dT~e~^=rf T . 
Jto h 

Interestingly, for 7 > 4/3, the observers following the fluid congruence will reach infinite ex- 
pansion in finite proper time. This is an indication that the fluid experiences some kind of 
singularity in the future [9] . 



4 Conclusion 



As an application of centre manifold theory, we have studied the centre manifold of the tilted 
Bianchi type III model. For the cases 7 = 1 and l < 7 < 2 this centre manifold is 3- and 
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2-dimensional, respectively. We analysed the late-time evolution and determined the decay 
rates for the variables (eg. ([25]) for 7 = 1, and eq. (|29p for 1 < 7 < 2). This analysis is a 
prime example of dynamical systems theory (in particular, centre manifold theory) appearing 
in applied mathematics and physics. 
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